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Abstract 

The nonperturbative Hamiltonian eigenvalue problem for bound states of a quantum field theory 
is formulated in terms of Dirac's light-front coordinates and then approximated by the exponential- 
ly operator technique of the many-body coupled-cluster method. This approximation eliminates any 
need for the usual approximation of Fock-space truncation. Instead, the exponentiated operator is 
truncated, and the terms retained are determined by a set of nonlinear integral equations. These 
equations are solved simultaneously with an effective eigenvalue problem in the valence sector, 
where the number of constituents is small. Matrix elements can be calculated, with extensions of 
techniques from standard coupled-cluster theory, to obtain form factors and other observables. 
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I. INTRODUCTION 



We wish to construct a method to compute hadron structure in terms of wave functions by 
solving the light-front Hamiltonian eigenvalue problem for quantum chromodynamics. The 
use of light-front quantization [l[ allows for a well-defined Fock expansion (2J, as indicated 
schematically in Fig. [T] for a proton. The fundamental eigenvalue problem 

(K.E. + V QC d)\p) = E p \p), (1.1) 

with K.E. representing the kinetic energy, Vqcd the interactions of QCD and E p = 
y/ m p + P 2 , is equivalent to a coupled set of integral equations, depicted in Fig. [2l when 
projected onto Fock sectors. 

In light-front quantization, the mass eigenvalue problem is 

V-\P)= p + M R), V\P) = P\P). (1.2) 

Here we define the light-front energy as P — [E — P z ), conjugate to the light-front time 
x + = t + z; and the light-front momentum P = (P + = E + P z , P± = (P x , P y )), conjugate 
to the spatial coordinates x = (x~ =t — z,x,y). Because p + is positive for all constituents, 
there are no spurious vacuum contributions to eigenstates; particles cannot be produced from 
the vacuum and still conserve p + . We also have a boost-invariant separation of internal and 
external momenta. Wave functions are then functions of longitudinal momentum fractions 
Xi = p\ /P + and relative transverse momenta ki± = pi± — XiP± of the constituents. 

To have a finite set of equations for the wave functions, there must be a truncation of some 
sort. The standard truncation is one of Fock space, to limit the number of constituents in- 
cluded in the calculation. Unfortunately, this introduces various complications, particularly 
uncanceled divergences [3| • For example, the Ward identity of gauge theories is destroyed by 
the truncation because a truncation of the number of photons in flight can eliminate vertex 
loops and self-energy loops on one leg but not the other Q]. We avoid these difficulties by 




FIG. 1. Fock-state expansion for a proton in terms of quarks and gluons. 




FIG. 2. Coupled equations for the Fock-state wave functions of a proton. 
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FIG. 3. The T operator. Each diagram corresponds to a contribution to T with emission of one or 
more bosons, represented by dashed lines, from the heavy fermion. The circled vertex is labeled 
by a function t n that determines the allocation of momentum to each boson. 

making a different sort of truncation, within the context of the light-front coupled-cluster 
(LFCC) method 0-0]. 

The LFCC method writes the eigenstate in the form \ip) = \J~Ze T \(j)) ) with \<f>) a valence 
state with the smallest number of constituents consistent with the desired quantum numbers, 
T an operator that adds constituents, and \J~Z a normalization factor. The operator T 
conserves J z , light-front momentum P, charge and any other quantum number that should 
be conserved. We then define an effective Hamiltonian V~ = e~ T V~e T and a projection P v 
onto the valence sector, and the original eigenvalue problem ( 11.2jl is equivalent to 

M 2 + P 2 

PvV~\<t>)= p + (1.3) 

and 

(l-P v )V=\<f>) = 0. (1.4) 

The second equation is an auxiliary equation for the operator T. 

At this point, the formulation is still exact and still reduces to an infinite set of equations 
for the infinite set of terms possible for T. The LFCC approximation is a restriction of T to 
one or a few terms and a corresponding restriction of 1 — P v to generate only enough auxiliary 
equations to solve for the terms kept in T. Many of the mathematical manipulations needed 
to carry out this approach can be found in the original many-body coupled-cluster method 0, 



101 ] . There the valence state is one of many particles, not a few, and the T operator introduces 
correlated excitations of the single-particle states without changing particle number. Despite 
the large difference in the physics, the mathematics is much the same. 



II. SAMPLE APPLICATION 

To illustrate the use of the LFCC method and to compare it with ordinary Fock-space 
truncation, we consider a light-front analog of the Greenberg-Schweber model [ll[]. I n the 
model, a heavy fermionic source emits and absorbs bosons without changing its spin. The 
full T operator would then be an expansion in multiple emissions of bosons, as shown in 
Fig. [3j The action of the exponentiation of T is given in Fig. HJ This is to be compared with 
a Fock-state expansion of the eigenstate, as presented in Fig. [5j where the wave function i[)q 
is equivalent to the LFCC valence-state amplitude 0. 

A truncation of the T operator to only the first term, single-boson emission, leaves only 
the t\ terms of Fig. |H The valence eigenvalue problem for and the auxiliary equation for 
t\ are then the coupled system of nonlinear equations depicted in Fig. EJ These are to be 
compared with the equation for the Fock wave function ip\ in Fig. [7] obtained by truncating 
the Fock-state expansion at two bosons and then integrating out the zero-boson and two- 
boson wave functions. A key distinction between the two approaches is in the self-energy 
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FIG. 4. The action of e T on the valence state \cj)). Factorials and other combinatoric factors are 
not shown. 
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FIG. 5. Fock-state expansion for the eigenstate The circles are labeled by the relevant wave 
function ip n . 



loop for the fermion; in the LFCC approach, it is the same everywhere it appears in Fig. |6j 
For the Fock-space truncation approach, it is different in each Fock sector, and the form 
that appears in Fig. [7] is dependent on the momentum of the spectator boson. 

In jHJ we show that the simplest truncation of T leads to an analytically solvable auxiliary 
equation and provides for this simple model the exact answer! The solution can be used to 
compute matrix elements, in particular the Dirac form factor for the fermion. Q 
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FIG. 6. Valence-sector eigenvalue problem and auxiliary equation in the LFCC method. The crosses 
represent kinetic energy contributions. The inhomogeneous term on the right of the auxiliary 
equation is the bare emission vertex from the original model Hamiltonian. 



1 The limitation of the model to a heavy fermion that cannot flip its spin eliminates consideration of the 
Pauli form factor; this is not a limitation of the LFCC method. 
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FIG. 7. Equation for the one-boson/one-fermion wave function ipi obtained in a Fock-space trun- 
cation to no more than two bosons. 



III. SUMMARY 



The LFCC method provides a nonperturbative approach to the solution of Hamiltonian 
eigenvalue problems in quantum field theories. It assumes that the Hamiltonian is suit- 
ably regulated, so that the integrals involved are rendered finite. The advantages of the 
method are that it does not require a Fock-space truncation and that it does not generate 
sector-dependent or spectator-dependent self-energies. Also, as an approximation, it is sys- 
tematically improvable through the addition of terms to the T operator, each constructed 
to add more particles to the valence state. 

An application to a simple model has been discussed here and more fully in [5|. An 
application to the dressed-electron state in QED is investigated in [gj, and additional work 
there is anticipated (p, 0]. The method should also be applicable to QCD, with light-front 
holographic QCD [12j as a natural starting point; there light-front holography may provide 
a useful approximation to the LFCC valence eigenvalue problem for QCD. 
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